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ABSTRACT 

f-f.' 

Thla,paper<-addresaes  the  problem  of  passive 
underwater  contact  tracking  In  the  presence  of 
mlsoodellng.  A  configuration  Is  proposed  which 
Includes  contact  state  estisMtion  and  hypothesis 
testing  coupled  with  an  evidential  reasoning 
capability  to  resolve  uncertainty  In  audellng 
assusiptlona.  The  Dempster-Shafer  theory  of  evidence 
is  applied  to  direct  the  generation  of  evidence 
required  to  classify  modeling  anomalies  and  provide 
compensation.  Experimental  results  are  presented 
that  desunstrate  the  reduction  of  modeling 
uncertainty  with  the  accumulation  of  evidence. 


I .  INTRODUCTION 

Contact  tracking  In  the  underwater  environment 
deals  prlsiarlly  with  processing  passively  observed 
measureoient  data  under  a  variety  of  operating 
conditions  [1].  The  efficient  extraction  of 
contact-related  parameters  from  measurement  data  Is 
contingent  upon  judicious  modeling  of  the  tracking 
envlronsient .  The  tracking  ondel  should  provide  an 
accurate  but  concise  description  of  the  underlying 
principles  and  assumptions  associated  with  contact, 
observer,  and  sensor  dynamics;  acoustic  sound 
propagation;  and  the  measurement  generation 
process.  Encumbrances  arise  when  abnormalities 
occur  within  the  estimation  process.  These 
abnormalities  are  exhibited  as  Irreconcilable  trends 
In  messurement  residuals,  indicative  of  mlsmodellng 
chat  has  occurred. 

Recent  work  has  addressed  the  detection  of 
modeling  anosiallos  through  analysis  of  measurement 
residuals.  In  particular,  hypothesis  testing 
techniques  have  been  employed  to  detect  bias  within 
temporally  delineated  segments  of  residuals  [2]. 
These  testing  techniques  have  been  Integrated  with 
contextual  Information  processing  In  s  prototype 
knowledge-based  system  for  targeting  solution 
Integration  and  assessment  (TARSIA).  TARSIA  [3] 
focused  on  the  intelligent  control  of  a  generalized 
maximum  likelihood  estimator  and  the  use  of 
numerically  derived  solution  attributes  as  the  basis 
for  heuristic  analysis  In  assessing  solution 
qtiallty.  As  such,  TARSIA  relies  on  synergistic 
coupling  between  numeric  and  symbolic  components  as 
Indicated  in  figure  1. 


Figure  1.  Functional  Components  of  TARSIA 


This  paper  extends  these  previous  efforts  by 
addressing  the  problem  of  uncertainty  In  the 
classification  of  data  and  modeling  anomalies  and 
Its  liapact  on  solution  quality  assessment  through 
the  application  of  evidential  reasoning  methods. 

Evidence  pertaining  to  the  nature  of  data  and 
modeling  anomalies  may  be  educed  from  several 
sources.  Including  Intelligence  and  cracking  data, 
sonar  Information,  ant^^sound  channel 
characterization.  Tc'  Is  amnifesCed  as  patterns 
within  time  series  (e.g.,  measurements  and 
measurestent  residuals)  and  trends  In  solutions. 

Often,  more  than  one  anomaly  type  might  be 
considered  based  on  a  given  sec  of  evidence.  This 
results  In  uncertainty  as  to  the  nature  of  the 
anomaly  and  a  need  for  additional  Information  to 
resolve  the  ambiguity. 

The  reasoning  technique  considered  for  Che 
interpretation  and  combination  of  evidence  la 
Dempster-Shafer  theory  [L,3].  An  Important  feature  3 
of  the  theory  la  Its  ability  to  narrow  the  set  of  Q 
likely  hypotheses  with  the  accumulation  of  evidence, 
while  retaining  the  notion  of  Ignorance  when  — 

insufficient  evidence  Is  presented.  This  Is  In 
contrast  with  Bayesian  techniques,  which  fall  to 
discriminate  between  equal  uncertainty  and  lack  of 
knowledge.  With  the  utilization  of  Dempster-Shafer 
theory,  a  set  of  plausible  anomaly  types  can  be  — 
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P(Ai)  -  1  -  S(Af)  . 


(4) 


foniiul«ted  that  allows  for  the  simultaneous  pursuit 
of  confirming  evidence  to  further  reduce  the  size  of 
the  set,  and  the  inclusion  of  candidate  anomaly 
types  in  the  formulation  of  the  parameter 
uncertainty  descriptions. 

A  brief  description  of  Dempster-Shafer  theory  is 
provided  in  the  next  section.  A  conceptual  model  of 
a  processing  configuration  that  incorporates 
Dempster-Shafer  theory  is  then  presented.  This  is 
followed  by  a  oiathematical  formulation  of  evidence 
generated  subsequent  to  data  processing. 

Experimental  results  based  on  angle-of-arrlval 
measurements  are  presented. 

II.  DEMPSTEE-SHAFER  THEORY 

The  Dempster-Shafer  theory  of  evidence  has  as  its 
basis  the  formulation  of  a  frame  of  dlscemsMnt. 

This  is  a  mutually  exclusive  and  exhaustive  set  of 
allowable  hypotheses,  denoted  by 

0  ■  {ai,  a2i  • ■ • ,  aL}  ,  (1) 

where  a^,  lwl,L  denote  the  individual  (singleton) 
hypotheses.  The  power  set  associated  with  6  is 
the  set  of  all  possible  subsets  of  6.  This  is 
illustrated  in  figure  2  for  a  four-eleoMnt  frasM  of 
discernment. 


(•«)  (stl  (*11  (•«! 


Figure  2.  The  Power  Set 

As  evidence  regarding  the  validity  of  certain 
hypotheses  becomes  available,  probability  mass  (or 
belief)  is  distributed  to  appropriate  subsets  within 
the  power  set  with  the  stipulation  that 

tB(Ai)  w  1  ,  (2) 

where  A^  represents  an  element  of  the  power  set, 
and  a(A^)  represents  the  non-negative  basic 
probability  assigtusents  (BPA).  Belief  not  assigned 
to  particular  subsets  is  assigned  to  8  and 
represents  uncertainty  in  the  interpretation  of 
evidence.  Associated  with  each  Aj^  is  an 
evidential  interval  [S,P]  where  S  represents  the 
support  for  the  element  and  P  represents  its 
plausibility  [6].  This  interval  corresponds  to  a 
minimum  and  maximum  likelihood  for  the  given 
subset.  When  evidence  confirms  a  subset,  S  is 
calculated  for  each  Aj  by  susaing  the  BPA  of  A{ 
and  all  its  subsets,  i.e., 

S(Ai)  -  2;  n(Aj)  .  (3) 

AjCAi 

Plausibility,  an  indication  of  the  lack  of 
disconf irming  evidence,  is  computed  as 


where  Af  denotes  the  complement  of  A^. 

Ignorance  regarding  the  validity  of  a  particular 
hypothesis  is  represented  by  the  evidential  interval 
[0,1],  while  [1,1]  and  [0,0]  represent  the  Boolean 
notions  of  true  and  false,  respectively. 

As  multiple  bodies  of  evidence  become  available, 
the  need  exists  to  merge  the  evidence  into  a  single 
belief  representation.  This  is  accomplished  through 
the  use  of  Dempster's  rule  of  combination,  given  by 

mi@m2(A£)  wK  ^  mi(X)m2(Y),  (5) 

xm-Ai 

K-1  -  1  -  F  m],(X)a2(Y)  .  (6) 

xnY-# 

Dempster's  combination  rule  is  both  cosmutative  and 
associative,  providing  an  order-independent  method 
for  merging  evidence.  Normalization  by  K  disallows 
any  probability  mass  being  assigned  to  the  null  set, 
since  this  would  violate  Che  exhaustive  set 
requirement  for  the  frame  of  discernment. 

A  prims  concern  when  applying  De^ster-Shafer 
theory  is  the  potential  computational  explosion  that 
occurs  when  the  else  of  the  frame  of  discernment 
becomes  large.  One  recent  approach  to  this  problem 
is  to  construct  multiple  frames  of  discernment 
linked  together  through  compatibility  relations  [7]. 
That  is,  in  lieu  of  defining  a  single  frame  that 
encompasses  all  possible  hypotheses,  the  hypothesis 
space  is  delimited  into  natural  partitions,  each  of 
which  can  be  addressed  by  a  reduced  order  frame.  A 
network  of  compatibility  mappings  is  then 
constructed  such  chat  evidence  producing  SPAs  in  one 
frame,  8j^,  results  in  analogous  BPAs  assigned  in 
another  frame,  8g.  A  compatibility  relation, 

8a, B>  defines  which  elements  in  two  distinct 
frames  are  true  sisiultaneoualy,  with 

®A,B  9  ®A  *  ®B  • 

Compatibility  mappings  are  given  by 

^A-^B^^k^  ”  {'>j|(*itbj)  C  8y^^B,  ai  C  A^}  ,  (8) 

where  at  least  one  pair  (a£,bj)  is  specified  for 
each  A^  in  8/^.  In  this  way,  any  frasie  may  be 
linked  to  any  other  frame  in  the  network. 

III.  PROPOSED  PROCESSING  CONFIGURATION 

A  configuration  for  integrating  Dempster-Shafer 
evidential  reasoning  with  statistical  parasieter 
estisiatioo  for  anosuly  classification  is  illustrated 
in  figure  3.  When  anosialies  are  detected,  the  data 
are  reprocessed  with  suspect  data  segamnts  omitted 
to  preclude  their  influence  on  the  generated 
solution.  They  are  oxinitored,  however,  for  clues  as 
to  the  nature  of  the  anomaly.  In  this  way,  residual 
characteristics  are  based  on  a  pivotal  unbiased 
solution  derived  from  a  consistent  data  set. 
Preliminary  evidence  is  generated  via 
characterisation  of  suspect  data  segments  in  terms 
of  their  perceived  features.  Refinesieot  of  these 
features  is  achieved  Iteratively  by  resegmenting 


length.  Further,  let  each  individual  measuremenr , 
Zj(k),  identified  by  measurement  type  j  and  time 
index  k,  be  governed  by  plant  and  measurement  models: 

x(k>ll)  ■  g[x(k),Y(k),k,k+a) ,  (10) 

*j(k)  »  f j[x(k',Y(k)l  +  nj(k),  (11) 

where  x(k)  la  an  n-dimensional  state  vector,  yCk) 
is  s  system  parameter  vector,  and  nj(k)  denotes 
a  sero-mean,  independent,  white  noise  sample  with 
N[0,o?].  The  MLS  yields  an  estimate  x  that 
minimizes  the  performance  index  (PI): 

PI  -  (1/N)(ZJ  D-1  Z^)  ,  (12) 


suspect  data  to  resolve  events  and  event 
boundaries.  Additional  evidence  is  furnished  by 
examining  solution  features.  These  include 
comparison  of  parameter  values  with  s  priori 
information  and  the  formulation  of  distance  measures 
between  alternative  solutions.  A  more  detailed 
description  of  these  sources  of  evidence  is  provided 
in  the  following  section. 

The  evidential  reasoning  component  contains  a 
natwork  of  frames  of  discernment,  where  a  frame  is 
constructed  for  each  measurement  type  under 
consideration.  Measurement  frames  contain  only 
those  anoawly  types  to  which  the  meaaureaient  in 
question  Is  sensitive.  A  hierarchically  related 
solution  frame  la  linked  to  each  measurement  frame; 
and,  as  evidence  is  applied  to  measurement  frames, 
compatibility  awpplngs  translate  the  evidence  to  the 
solution  frame.  In  this  way,  evidence  across 
measurement  sets  is  merged.  Solution  extracted 
evidence  is  applied  to  the  solution  frame. 

Inferences  regarding  anomaly  classification  are 
produced  primarily  through  examination  of  the 
solution  fraaie. 

The  classification  of  anomalies  generally  results 
in  modification  of  modeling  assumptions  over 
segments  of  data.  In  the  event  that  the  anomaly 
type  remains  unresolved,  a  plausibility  ranking 
motivates  the  search  for  confirming  evidence  or,  if 
this  is  unattainable,  directs  the  inclusion  of  these 
anooialy  types  into  the  formulation  of  uncertainty 
region  descriptions. 

IV.  MATHEMATICAL  FORMULATION  OF  EVIDENCE 

The  estimator  represented  in  figure  1,  an 
iterative  maximum  likelihood  estimator  (HLE)  used 
for  mapping  data  Into  solutions  and  integrating 
solutions,  is  described  by  Murphy  [8].  Its  outputs 
include  the  solution  In  the  specified  coordinate 
system  and  a  square  root  realisation  of  the  Fisher 
information  matrix  (FIM).  Based  on  these,  various 
pieces  of  evidence  are  created.  To  facilitate 
subsequent  discussion,  let  Z  be  defined  as  an 
N-dimensional  measurement  vector: 


where  Z^  denotes  the  measurement  residual  vector 
with  elements 

z^jCk)  »  Zj(k)  -  f j(x(k),Y(k)l  ,  (13) 

and 


D  ■  E{njn5}  ■  Diag{iJ^}.  (14) 

Note  that  for  unbiased  estimates  the  PI  becomes 
X^-dlstributed  with  an  expected  value  of  unity. 

In  a  similar  fashion,  the  data  partitions  of 
equation  (9)  yield  local  values  of  the  PI. 


Solution  Consistency 

The  priaiary  attributes  used  to  ascertain  the 
conaistency  of  solutions  include  both  local  and 
global  assessments  of  the  PI,  and  the  output  of  a 
hypothesis  testing  technique  in  which  messurement 
residuals  are  interrogated  for  bias.  In  the  former 
case,  thresholds  are  set  to  detect  the  occurrence  of 
improbable  values  of  PI  based  on  a  specified 
probability  of  false  alarm.  In  the  latter  cate,  two 
hypotheses  are  formulated.  The  null  hypothesis  Hq 
assumes  that  the  residuals  within  a  given  data 
segment  are  zero  mean.  Under  the  alternative 
hypothesis  H^,  the  data  residuals  exhibit  a 
significant  bias.  That  is, 

Hq:  “  (9fj/9x(k))  5x(k)  ♦  nj(k)  , 

Hp  z,j(k)  »  [3fj/3x(k)}  5x(k) 

♦  aj(k)  ♦  nj(k), 

where 

5x(k)  -  x(k)  -  x(k)  , 

and  mj(k)  denotes  a  least-squares  fit  to  the 
residuals  within  the  data  segment.  The 
Neyman-Pcarson  test  procedure  is  applied  [2]  to 
yield  the  test  criteria 


(15) 


(16) 


(in 


where  Z^,  lwl,L,  denote  partitions  of  srbitrary 


where  X  is  chosen  to  satisfy  a  specified 
probability  of  false  alarm.  A  data  segment  is 
deesied  consistent  when  the  hypothesis  test  is 
passed,  indicating  unbiased  residuals  in  accordance 
with  the  theoretical  model.  In  general,  bias 
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exhibited  within  one  or  more  data  segments  can  cauae 
the  solution  x  to  becosie  biased  and.  in  the  process 
of  minimization,  results  in  other  data  segments 
appearing  biased.  Thus,  when  a  particular  data 
segment  la  deemed  to  be  the  source  of  modeling 
inconsistency,  it  is  removed  from  the  processed  data 
set,  while  its  associated  attributes  continue  to  be 
monitored.  In  this  way,  features  extracted  from 
monitored  segments  are  baaed  on  consistent  solutions 
and  provide  evidence  as  to  the  nature  of  detected 
anomalies. 


is  the  state  transition  matrix.  Note  that  when  x  is 
unbiased,  the  expected  value  of  the  residual  becomes 
that  of  the  modeling  error  mapped  into  measurement 
space.  Depending  on  operating  conditions,  the 
modeling  error  can  easily  dominate  the  residuals. 
Hence,  for  a  given  class  of  modeling  errors,  the 
effect  on  measureiaent  residuals  is  predictable. 

Note,  however,  that  modeling  parameters  are  subject 
to  observability  requireoients  just  as  are  state 
parameters. 


Perturbation  Analysis 

Given  an  estisiate  x  deemed  consistent  prior  to 
the  detection  of  an  anomaly,  the  effects  of 
mismodellng  on  the  residuals  is  now  examined. 
Equations  (10)  and  (11)  can  be  approximated  in  a 
Taylor  series  expansion  as 

x(k>l)  »  g[x(k),  Y(k)t  k,  k+U 

+  (3g/3i(k)l  Sx(k)  ♦  (3g/3Y(k)l  SY(k)  .  (18) 
ij(k)  a  fj[x(k),  Y(k)l  ♦  [3fj/3x(k)]  8x(k) 

♦  (3fj/3Y(k)l  «Y(k)  ,  (19) 

where 

JY(k)  ■  Y(k)  -  Y<k)  .  (20) 

This  yields  the  approximate  residual  sequence 
**j(k)  a  nj(k)  ♦  [3fj/3x(k)]  4x(k) 


Solution  Features 

Solutions  derived  from  independent  data  sets  can 
be  integrated  by  forming  and  processing  pseudo¬ 
measurement  and  pseudo-residual  sets  of  the  form 

^i  ■  >1  *i  •  (25) 

Zei  -  Ritxi  -fi(i)l  ,  (26) 

where  xj^  and  denote  local  solutions  and  their 
square  root  information  matrices,  respectively,  and 
X  denotes  the  global  (integrated)  solution,  x^  and 
X  need  not  be  of  the  same  dimension  or  coordinate 
system  as  indicated  by  the  swpplng  function  f^. 

Note  that  equation  (26)  provides  a  measure  for 
evaluating  solution  cluatering  as  formulated  by 
Llndgran  at  al.  [9].  Threshold  comparisons  of  the 
local  PI  in  this  case  are  used  to  provide  evidence 
regarding  the  similarity  of  candidate  solutions  when 
the  individual  solutions  are  based  on  independent 
sndeling  assumptions. 


♦  (3fj/3Y(k)l  SY(k)  .  (21) 

Consider  the  case  where  system  parasMter 
mismodellng  comsiences  at  time  kg.  The  effect  of 
modeling  changes  after  time  kg  can  be  computed  from 


3f. 


z_,(k_+l)  »  n,(k_+l)  +  — 


ej'  m 


j 


3x(k>l) 

m 


♦[k_,  k  +1.  Y(k_),  x(k_))  6x(k_) 


R  3f , 
t  E  ^ 


i-1  3x(k  ♦!) 


♦(k^,  k^+i.  Y(k,*i).  x(k^*l)) 


3Y(k^>i) 


6Y(k  ♦!). 


(22) 


Further  evidence  can  be  generated  through  a 
partitioning  of  the  pseudo-measurements  into  reduced 
dimension  components  and  forming  partitioned 
distance  measures  of  the  form 


^ei  “  [*il^*ir^il^*^^  :  "i2^*i2"‘^12^*^^  ;  •'•  ]■ 

(27) 

Since  is  non-diagonal,  this  requires  formulation 
of  marginal  probability  density  descriptions  that 
can  be  mapped  into  reduced-order  square  root 
information  matrices. 

Examples;  Antle-of-Arrival  MeasuremeT.i.s 

The  effects  of  certain  anoeialy  types  on 
angle-of-arrival  measurements  are  now  provided. 
Consider  the  measurement  types  bearing  and 
depression/elevation  (D/E)  a.-tgle  given  by 


where 

6x(k4'l)  ■  [3g/3x(k)I  6x(k) 

♦  (3g/3Y(k))  6Y(k)  ,  (23) 

and 

♦Ik.  kU,  Y(k),  i(k)]  -  kt 

3x(k) 


B  ■  tan"E(rg/ry)  ,  (28) 

Sp/g  ■  tan“E(r,//rJ»r^)  ,  (29) 

where  x(k)  is  de'.ned  in  a  Cartesian  coordinate 
system  with 

XT  “  lrxt»  fyt»  ’'xt*  ''yt*  ''ztl^  • 

*o  "  t^xo*  '^yo»  '^xo*  '^yo*  *^zo^^  • 
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X  .  [fj,,  Ty,  V,,  Vy,  r,]T  ,  XT  -  X(,  , 

'l  0  T  0  0‘ 

0  I  0  T  0 

♦(k,k+l)  «  0  0  I  0  0  .  (30) 

0  0  0  1  0 

0  0  0  0  1 

and  where  X(,  ilg,  and  S  denote  target,  obaerver, 
and  relative  state  vectors,  respectively.  Here, 
r,(  Is  defined  as  the  depth  of  a  target  image  when 
reflections  of  sound  occurs  on  the  bottom  or  surface 
of  the  sound  channel.  Also  consider  the  conical 
angle  sieaaurement  type  as  observed  with  a  linear 
array: 

0^  -  tan“^r|,/-/r?Tr|),  (31) 

where  r,|  and  denote  parallel  and 
perpendicvlar  displacements,  respectively,  of  the 
target  from  the  linear  array.  The  modeling 
perturbation  Syfk)  is  now  defined  for  certain 
anomalies.  Consider  first  an  unsndeled  target 
maneuver.  Often,  a  constant-velocity  assumption  is 
placed  on  target  motion,  and  a  change  in  target 
course  or  speed  violates  this  assumption.  This  can 
be  related  as  plant  inputs 

v^((k4'l)  «  Vx((k)  >  ,  (32) 

Vy((k>l)  ■  Vy((k)  ♦  Avy  ,  (33) 

where  and  dvy  are  normally  aero  and,  for 
an  instantaneous  change  in  velocity,  are  non-sero 
only  at  time  k^,.  Similarly,  consider  an  unmodeled 
change  in  relative  image  depth  dr,  occurring  at 
tisw  kg  due  to  either  varying  bottom  conditions,  a 
change  in  propagation  path,  or  a  target  depth 
change.  These  modeling  perturbations  can  be 
represented  as 

Sylk)  ■  (dv,  dvy  dr,lT  (for  k  -  kg)  .  (34) 

Finally,  consider  an  additive  bias  placed  on  the 
respective  measurement  types  and  cooencing  at  time 
kg.  Here 

8Y(k)  -  dsj(k)  (for  k  >  kg).  (35) 

Then,  it  is  a  straightforward  calculation  to  show 
that 

dS(k)  -  daQ(k)  *■  [dB(k)](k  -  kg)T  , 
d0D/E('‘^  ■  dX5/g(k) 

♦  dr,(k)  cos  dD/g(k)/r,(k)  -  [d0[)/E(k))(k  -  kg)T, 
de(;(k)  •  dsc(k) 

♦  dr,(k)  sin  0Q/£(k)/(r,(k)  tan 

♦  [d^(k))(k  -  kg)T  ,  (36) 

where 

dA(k)  ■  [dvg  cos  a(k)  -  dVy  sin  B(k)]/r(k)  , 


d6Q/£(k)  *  [dv^  sin  B(k)  *  dVy  cos  B(k)] 

•  sin  • 

de^fk)  *  [drj^  cos  ^(k)  cos  6c(k) 

-  dr,!  sin  ec(k)J/r,(k)  , 
r(k)  -  (r2(k)  ♦  r2(k)ll/2  , 
r,(k)  -  [r2(k)  ♦  r2(k)ll/2  . 

?(k)  »  tan"'^(r,(k)/r^(k) J  ,  (37) 

and  where  solution  error  and  measurement  noise  are 
asausmd  to  be  dominated  by  the  modeling  error. 

Note  that  the  perturbations  consist  of 
measurement  bias  and  geometry-dependent  components. 
Furthermore,  image  depth  change  is  constant  after 
time  kg,  yet  its  effect  on  measurements  is 
geometry-dependent.  The  target  maneuver,  however, 
while  also  modulated  by  the  geosietry,  has  a  linear 
variation  in  time.  These  features  can  be  postulated 
as  Jumps,  drifts,  and  curvature  in  the  measurement 
residuals. 

If  independent  solutions  are  produced  based  on 
data  sets  before  and  after  time  kg,  it  is  seen  that 

Xa(kg)  -  Xp(kg)  •  to  0  dVg  dVy  dt , )  ^  (38) 

for  the  target  swneuver  and  image  depth  change.  The 
effect  of  bias  is  unpredictable  and  depends  on  the 
precise  nature  of  the  bias.  Bias  removal  or 
compensation  techniques  can  be  evaluated  by 
examining  how  well  solutions  cluster  with 
conipensation.  Finally,  if  image  depth  is  not 
estimated  but  presuawd  known,  then  modeling  error 
will  contaminate  the  other  target  parameters  and  the 
effect  is  again  unpredictable.  Once  nx>re,  solution 
clustering  can  be  assessed  to  validate  the  assumed 
image  depth 


V.  EXPERIMENTAL  RESULTS 

Experimental  results  are  presented  to  demonstrate 
the  utility  of  the  proposed  concept.  The  network  of 
frasws  shown  in  figure  4  was  constructed  with  the 
frasws  of  discernment  defined  in  tables  1  and  2. 
Static  bias  is  defined  as  occurring  instantaneously 
and  maintained  constant  over  a  segment  while 
dynamic  bias  is  defined  as  tisie-varying  but  with  no 
discontinuities  in  the  waveform.  Since  bias  can 
occur  in  any  Masuresmnt  set,  both  static  and 
djmamic  bias  for  each  measuresient  set  is  linked  to 
systematic  error  in  the  solution  frame.  Iisage  depth 
changes  are  presumed  to  occur  instantaneously. 

Sensor  alssx>deling  applies  only  to  the  linear  array 
and  solutions  derived  from  conical  angle  data. 

Basic  probability  assignments  were  developed  based 
on  the  anomaly  set  defined  above,  with  separate  BPAs 
defined  for  interpreting  the  features  of  measureswnt 
residuals  for  each  respective  measurement  type. 
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Figure  4.  Fraae  Hierarchy 


Table  1.  Anooaliea  Conaidered 


Hypothaaia _ Deacriotion 


SB 

DB 

TM 

BZ 

SM 

SB 

Static  Bias 

Dynaaic  Bias 

Target  Maneuver 

Image  Depth  Change 

Sensor  Nlsaodellng 
Systeaatic  Error 

Table  3.  Frai 

aes  of  Oiscammant 

Blamenta 

Compatibility  Relations 

Frame 

of  Fraaa 

with  Solution  Frame 

Solution 

{SI,Tn,BZ,SM} 

— 

Bearing 

{SB,DB,TM> 

(SB,S8),(0B,SB),(1M,TM) 

D/E  Angle 

{SB,DB,TM,RZ} 

(SB,SE),(DB,SE),(TM,'m), 

(RZ,RZ) 

Conical 

{SB,DB,'m,RZ,SN>  (SB,SE),(DB,SE),(TH,ITf) 

Angle 

(RZ,RZ),(SN,SN) 

Three  anoawly  typea  were  aiaulated  in  cheae 
experiflwnta:  target  aaneuver,  static  bias  in  0/E 
angle,  and  iaege  depth  change.  Test  procedures  to 
generate  confinalag  evidence  by  coapensating  for 
these  anoealy  typea  were  developed  with  subsequent 
BFA  distribution.  The  geoaatry  considered  is  shown 
in  figure  5,  where  the  naneuvering  target  trajectory 
Is  Indicated  by  a  daabed  line.  Exaaples  of  residual 
features  in  angle-of-errival  ■eesurewnt  residuals 
that  are  characteristic  of  these  foms  of 
oiseodeling  are  shown  in  figitre  6.  Experiiaeotal 
results  for  the  case  of  the  target  aaneuver  are 
shown  in  figure  7,  where  support  for  that  anooaly 
type  is  plotted  as  a  function  of  the  site  of  the 
data  window  following  anonaly  detection  that  is 
available  for  feature  extraction.  Support  is 
plotted  for  low,  nediusi,  and  high  noise  cases. 

Figure  7a  depicts  the  belief  in  target  aaneuver  as 
the  exhibited  anoaaly  type  when  only  bearing  is 
processed,  while  figure  7b  depicts  the  case  where 


bearing  and  D/E  angle  are  processed  jointly.  The 
target  maneuver  presents  drift  characteristics  that 
become  increasingly  apparent  with  time.  This  is 
reflected  in  figure  7  in  the  convergence  of  support 
for  this  anoaaly  type. 

Figure  8  describes  results  obtained  when  static 
bias  in  D/E  angle  and  iaage  depth  change  are 
considered.  In  these  cases,  the  anosialy  types 
exhibit  nearly  identical  features  at  the  occurrence 
of  the  anoaaly  and  become  only  marginally 
discemable,  given  sufficient  time  and  a  suitable 
tracking  geometry.  Hence,  figure  8  depicts  support 
for  the  individual  anomaly  types  as  well  as  their 
union.  Figure  8a  describes  results  when  static  bias 
is  the  actual  anoaaly,  and  figure  8b  describes  the 
image  depth  change  case.  These  results  accentuate 
the  utility  of  Deapster-Shafer  theory  for 
representing  modeling  uncertainty. 


VI.  SUmUlT  AKD  CONCLUSIONS 

Iba  integration  of  the  Deapster-Shafer  theory  of 
evidence  with  statistical  parameter  estimation 
techniques  has  been  demonstrated.  Doing  this 
requires  development  of  a  f raae-of-dlscemaeot 
hierarchy,  feature  extraction,  basic  probability 
asaignaents,  and  formulation  of  methods  of  searching 
for  confiraing  evidence.  In  addition,  the  utility 
of  forming  a  plausible  set  of  anomaly  types  from 
wbicb  to  build  inferences  regarding  solution  quality 
has  been  demonstrated.  That  is,  compensation 
teebniquea  should  accoMOdata  all  anomaly  types 
deemed  plausible,  given  that  anomalies  are  detected. 

It  should  be  noted  that  while  BFA  distribution 
for  the  experiaents  presented  was  based  on  analysis 
of  evidence,  it  can  still  be  priswrily  subjective  in 
nature.  Deapster-Shafer  theory  accoMOdates  the  use 
of  known  a  priori  statistics,  but  also  allows  for 
conservative  distribution  of  belief  as  a  function  of 
the  integrity  of  the  evidence.  That  is,  by 
assigning  nore  likelihood  to  the  fraae  of 
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Figura  6.  Meaauraiwat  Parturbatlon  Exaaplaa 


1.  MMIM  OKT 


diacernaant  itaalf,  the  onua  for  hypothesis 
resolution  la  placed  on  the  evidence  gathering 
function.  Slallarly,  a  priori  Information  obtained 
In-sltu  should  be  allowed  to  effect  the  distribution 
of  belief. 

The  slaipllfied  examples  presented  are  Indicative 
of  the  processing  capabilities  required  to  generate 
and  Interpret  evidence  for  the  general  case  of 
parameter  estimation  In  the  presence  of  modeling 
uncertainty.  Feature  extraction  should  be  matched 
to  the  modeled  frame  of  dlscemsmnt  and  can  require 
a  high  degree  of  aophlstlratlon.  Temporal  reasoning 
capabilities  may  be  required  to  ascertain  event 
occurrence  and  event  boundaries.  This  is  especially 
true  when  multiple  events  can  occur  either 
almultaneoualy  or  sequentially.  Future  efforts  will 
address  this  latter  aspect. 
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